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ABSTRACT 


The use of atmospheric drag for slowing satellites in high 
energy, high apogee orbits to a lower energy, lower apogee orbit 
about the earth is investigated- The high energy orbit is assumed 
to intercept the earth's atmosphere- Guidance Tor the atmospheric 
phase of the trajectory may be done using the aerodynamic forces 
generated by the passage through the atmosphere- This research 
was concerned with the investigation of several methods of 
guidance during the atmospheric phase to cause a significant 
reduction in the final velocity as the vehicle leaves the 
atmosphere. In addition, the velocity direction was controlled to 
exit to a desired target orbit. Lastly excess aerodynamic lift 
was used to effect a plane change between the entry orbit plane 
and the exit orbit plane to achieve a desired orbit plane. 

The guidance methods were applied to a 3 degrees-of— freedom 
simulation which included an oblate earth gravity model and a 
rotating atmosphere. Simulation results were compared on the 
basis of speed of computation of the guidance parameters and 
amount of added velocity necessary to achieve the desired orbit. 
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Guidance For /n Aeroassisted Orbital Transfer Vehicle 


There is currently underway a study concerning the use of 
atmospheric braking for vehicle transfer from a high energy orbit, 
to a lower energy crbit about. the earth. The use of atmospheric 
drag is seen as an excellent way to decrease the vehicle speed in 
the high energy orbit to the speed required for the lower energy 
orbit. In addition to the decrease in the speed of the vehicle, 
the plane of the new orbit may not coincide with the plane of the 
original orbit. This study is concerned with the investigation of 
the guidance requirements for the atmospheric braking portion of 
the vehicle trajectory with targeting to a final orbit. 

Problems of this type are generally analyzed by the iterative 
methods of perturb«<tion analysis. The trouble with this case is 
that the time required for the calculation of the solution to the 
two point boundary Value problem arising from the perturbation 
analysis may be too time-consuming for use in the vehicle guidance 
system. The atmospheric flight phase of the trajectory is 
characterized by sensitivity to the time in which maneuvers are 
accomplished. This sensitivity may even lead to atmospheric 
capture of the vehicle or to a final orbit which is nowhere near 
the target orbit. In addition to the time sensitivity, the 
atmospheric flight is characterized by large changes in 

atmospheric densit}' and dynamic pressure on the vehicle. Due to 
the requirement that the guidance system operate the guidance loop 
rapidly and effectively, a reduced complexity set of differential 
equations was chosen for the guidance system model and approximate 
solutions were employed for guidance and targeting. These 
solutions were applied to a 3 degree of freedom, point mass 
simulation with full complexity simulation differential equations 
including an atmospheric model with perturbations which rotates 
with the earth and an oblate earth model with the J2 harmonic 
terms included in the gravity model. 

The figure below represents a shuttle launched experiment 
for testing the vehicle and the guidance system. The shuttle 
launches the vehicle into a a high apogee lob trajectory. Just 
before returning to the original launch point, the on— board engine 
increases the speed and changes the velocity direction to simulate 
a vehicle returning from a geosynchronous orbit. The new orbit 
intersects the "reasonable** atmosphere of the earth at an altitude 
of 400000 feet, passes through the atmospheric transfer portion of 
the trajectory during which atmospheric forces are employed to 
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effect velocity and orbital plane changes to target to the final 
desired orbit at the exit from the atmosphere. The target orbit 
is characterized by an apogee altitude of 200 miles and a specific 
orbital plane inclination. At the apogee of the achieved final 
orbit, additional changes are made in the orbit to circularize the 
final orbit and to reach the desired orbital inclination. For 
this simulation, the specified entry conditions and the desired 
exit conditions are tabulated below. 


ENTRY CONDITIONS 

H-400000 FT 
V-33828 FPS 
4.3 DEG 
i-27. DEG 


EXIT CONDITIONS 

H-400000 FT 
HP— 50. MILES 
HA— 200 . MILES 
i— 28 . 5 DEG 


FIGURE 1. PROBLEM DEFINITION IN GRAPHICAL FORM 


aerobrakinq ExPtniMtHt 


Station 100 (X - 100 in) 
Rake Plane 


\ 


cti 



/ 




ORIGINAL PAGE IS 
OF POOR QUALITY 



During the atmospheric transier portion ol' the trajectory, 
the AOTV is assumed to maintain a constant angle or attack with 
the resulting constant lift coelTicient CC L > and constant drag 
coefficient <C D >. Control of the AOTV during this phase of the 
trajectory is done by rolling the vehicle to orient the lift force 
to cause in-plane vertical acceleration for altitude control and 
out-of-plane transverse acceleration for control ol* the orbital 
plane inclination. The equations of motion used in the guidance 
model may be written as follows: 

GUIDANCE EQUATIONS OF MOTION 


d V 
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c D s 

p V| _ G d « 

dt 
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V dt 
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In the equations above, V is the inertial velocity magnitude, 
Vjg is the airspeed or speed relative to the rotating atmosphere, 
p is the density of the atmosphere at altitude H above the surface 
of the earth, G is the gravitational acceleration at altitude H 
above the surface, R is the distance from the center of the earth 
to the AOTV, M is the constant mass of the AOTV, C^ and C^ are the 
lift and drag coefficients of the AOTV, and <f> is ihR roll angle 
used to control the vehicle. Using the second equation, the 
guidance equation may be specified as follows: 

GUIDANCE ROLL COMMAND EQUATION 


C0SC4>> 


COM 


2 M 
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Two items are apparent from the equation above; Cl) The 
roll angle is specified in cosine form, and C2) The roll angle 
cosine depends upon the specification ol' the vertical 
acceleration. For the first item, use the guidance roll command 
equation to specify the cosine of the roll angle. Use the 
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direction of the desired plane change maneuver to determine the 
sign of the sine of the roll angle. The combination of the roll 
angle sine and cosine will then uniquely determine the roll angle. 
A characteristic of this method is that the roll angle "ratchets" 
in that the amount of roll angle for plane control may be too 
great and the sign of the sine may need to be reversed 
periodically to stay within an allowable deadband of the desired 
final inclination angle. Typical trajectories will display four 
or more sign reversals on the sine term during the atmospheric 
pass. The second item, the commanded vertical acceleration in H 
is a different story. This study explored four methods of 
completing the roll command equation. 

METHOD 1. ASSUME A CONSTANT ALTITUDE RATE 


For a constant altitude rate, the commanded acceleration is 
zero which simplifies the guidance equation. For the entry phase, 
this application is straightforward as the initial conditions to 
start the "equalibrium glide" are specified. For the exit phase, 
the initial conditions are not specified and one must use 
iteration to determine the speed with which to begin the exit 
phase and the altitude rate command for the exit phase, which in 
combination yield the desired exit conditions. CSee C. Cerimele’s 
paper listed in the references for a detailed explanation of this 
method. > This method does succeed but relies upon the iteration 
technique for commands during the exit phase. In actual practice, 
the commanded roll angle is augmented by commands which will 
produce desired values of altitude rate and dynamic pressure, as 
otherwise the command would drive the AOTV down the constant 
altitude rate path into the lower atmosphere from which it would 
not be able to initiate an exit phase. 


CONSTANT ALTITUDE RATE ROLL COMMAND 

COSCtf»> GOI1 “ 2 M IG - - GH ( 

C L S pV 2 R 


dH 

ar 


_dH 

ar COM 


GQ 1QB-0B com >] 


The GQ and GH terms are dynamic pressure and altitude rate gains 
used in combination with altitude rate, dH , and dynamic pressure 

ar 

, QB, to produce the desired trajectory upon entry. For the exit 

HH 

trajectory, GQ is zero and the is determined iteratively. 

dt 
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METHOD 2. ASSUME A CONSTANT ROLL ANGLE COSINE 

This is another easily implemented option as the constant 
roll angle cosine eliminates any control determination for the 
vertical acceleration. The sign of the roll angle sine is still 
chosen to produce the desired plane change. The determination of 
the actual roll angle cosine is done in an iterative manner using 
a numerical perturbation in the constant roll angle to determine 
numerical partials of the target conditions with respect to the 

chosen roll angle cosine. Again, this method will work, but it 

depends upon solving a two point boundary value problem 

iteratively, and speed of guidance calculations may be its demise. 

CONSTANT ROLL ANGLE COSINE COMMAND 

COSC^> COM » CONSTANT. 

METHOD 3. ASSUME AN ALTITUDE TIME HISTORY 

For this method, assume a polynomial in time for the altitude. 

Include in the polynomial enough coefficients to be determined to 

allow the method to fit an optimum criteria. For example: 

ASSUMED ALTITUDE POLYNOMIAL 

1. ASSUME H - H ± + d H ± T + d 2 H ± T 2 + d 3 H T 3 

ar 3T 2 I Ht 3 6 

d H ■ d Hj + d 2 Hj T + d 3 H T 2 
dt dt dt 2 d t 3 2 

d 3 H » 2 t d H. 3 - d H. - d 2 H. I’ 3 

at 3 r 2 ar ar ar 2 

t| [ 1 d 2 H 1 3 + T 2 t 2 d H ± + 1 d H.,3 + L H 2 3 = O 

Z 3t 2 3 dt 3 dt 

1’ 2 - - B - l B 2 - 4 A C 3 122 
2~S 

GUESS UNICNOWN VALUES AND ITERATE TO A SOLUTION. 
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d 2 H 
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The method above has merit, in that, it, will allow targeting to a 
specific altitude with time, altitude rate, and vertical 
acceleration involved in the iteration. Analytical approximations 
allow quadrature approximations Tor integration which allow speed 
increases in iterative passes. It is still an iterative scheme 
but is a fast iterative scheme with explicit targeting. 

METHOD 4. ASSUME A VERTICAL ACCELERATION FUNCTION 

This last method involves assuming that the AOTV will act 
like a damped harmonic oscillator with specil'ied damping ratio and 
natural frequency. The vertical acceleration command is then 
easily determined from the following equation: 

ASSUMED OSCILLATOR FORM 


d2 “ COM " 2 d H h goM 5 

3U 2 3t 

In t*he oscillator equation above, the selection or damping ratio, 
and natural frequency, are equavalent to specifying the 

manner in which the AOTV system will react* The appearance of 
H C qji the equation allows targeting to a specific altitude on 
the entry phase and assures that the altitude rate will be nulled 
at that altitude* For the exit phase, this same Form can be used 
with different system characteristics* 

COSC^ COM - ™ t C f!”> C 0M + G - — 2 J 

C L S pVg dt. 45 R 

The outstanding benefits of the guidance scheme listed above 
is that it is fast Cwith no iteration), it allows altitude 
targeting, and the same form of guidance can be used in the entry 
phase and the exit phase with appropriate changes in the system 
characteristics. 
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The figures shown below represent. the trajectory shape for the 
atmospheric phase, with the attendant dynamic pressure, and a roll 
cosine profile. The roll cosine does not exhibit the ’’ratcheting'’ 
which the roll angle does as the "ratcheting" is due to sign 
reversals on the sine function. 
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